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ABSTRACT 



We report a periodicity of ~1 day in the highly elevated X-ray emission 
from the protostar V1647 Ori during its two recent multiple-year outbursts of 
mass accretion. This periodicity is indicative of protostellar rotation at near- 
breakup speed. Modeling of the phased X-ray light curve indicates the high- 
temperature (~50 MK), X-ray-emitting plasma, which is most likely heated 
by accretion-induced magnetic reconnection, resides in dense (>5xl0 10 cm" 3 ), 
pancake-shaped magnetic footprints where the accretion stream feeds the new- 
born star. The sustained X-ray periodicity of V1647 Ori demonstrates that such 
protostellar magnetospheric accretion configurations can be stable over timescales 
of years. 
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Introduction 



Strong X-ray emission and collimated jets from newborn stars, so-called protostars, 
indicate that energetic magnetic activity plays an important role in star formation. However, 
the thick gaseous envelopes that obscure protostars at visual wavelengths complicate efforts 
to probe their innermost active regions. Low-mass stars like the Sun in the formation phase 
(hereafter, young stellar objects or YSOs) gradually accumulate mass from the parent cloud 
before igniting nuclear fusion in the stellar core. Matter falling from the cloud first forms 
a disk around the central star and matter in the innermost disk gradually accretes onto it. 
The young star's gravitational potential can accelerate infalling matter up to a few hundred 
km s _1 . When the accreting material collides with the stellar surface, it is shock-heated to 



a few MK; this th ermalized plasma can emit soft (E < 1 keV) X-rays (IKastner et al.l 12002 
Sacco et al.lboioh . 



Matter cannot fall from the disk onto the central star unless it loses a significant fraction 
of its angular momentum. Coupling of the large scale magnetic fields of the cent ral star to 



those of the innermost dis k is suspected to prompt the momentum transfer (e.g., IShu et al. 



19961 ; lHayashi et al.lll996l ). Magnetic reconnections triggered by the interactions may be 
involved in the ejection of a fraction of the infalling disk matter, along with most of the 
angular momentum, out of the system; the rest is accreted onto the central star. Bipolar jets 
or collimated outflows seen ubiquitously outside of the YSO's e nvelopes are likely laun ched 
by magnetohydrodynamic or magnetocentrifugal processes (e.g.. lReipurth fc BaHyll200ll . and 
references therein). 

Although protostellar magnetospheric accretion models are widely accepted, and have 
gaine d support from obser vations of specific, relatively evolved objects (e.g., classical T Tauri 
stars. iDonati et al.ll2011al Jbl). the accretion geometry - and hence the validity and applicabil- 
ity of such models at earlier protostellar evo lutionary stages - remai ns uncertain in the case 
of younger, more deeply embedded objects (IJohns-Krull et al.l 120091 ) . X-ray observations of 
rapidly accreting objects may offer a means to probe this magnetospheric accretion process. 
The bulk of a star's mass is accreted during very early (observationally, Class and I) phases, 
when the star is still deeply embedded in its parent cloud . Such protostars ten d to emit hard 
(>2 keV) X-rays with occasional rapid (<1 day) flares (jlmanishi et al.ll200lf ). Some of the 
flares have been proposed as arising in large magnetic loops that connect the inner accre- 
tion disk and the stellar surfa ce (i.e., in a non-solar- type magnetic structure), which may b e 
loaded by accreting material (ITsuboi et al.lll998t iMontmerle et al.ll2000t iFavata et al.ll2005l ). 
However, systematic studies of X-ray-emitting pre-main sequence stars in the Orion nebula 
suggest that the occurrence of the largest magnetic loop s, associated with ma jo r flare events, 
is not dependent on the presence of circumstellar disks ( Getman et al. 2008al Jbl; lAarnio et al. 
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Over the last decade, a small number of YSOs have been found to display sudden in- 
creas es in mass accretion rate, by as much as a few orders of magnitude (IHartmann fc Kenyon 
19961 ). Such eruptive YSOs are crudely classified as either FU Ori (FUor) or EX Lupi (EXor) 
types, wherein the former generally display major outbursts over timescales of decades, and 
the latter generally display smaller, shorter-duration, outbursts. One such eruptive YSO, 
V1647 O ri — which does not fit neatly into either the FUor or EXor eruptive class (see discus- 



sion m 



Teets et al.ll201ll ) — went into outburst during the period 2004-2006 and then, again 



in 20 08 (the latter eruption is evidently still ongoing). This Class I YSO (IMuzerolle et al. 



20051 . Principe et al. 2012, in prep.) offered the first direct evidence that X-ray activit y 



surges with increases in mass accretion activity (IKastner et al.ll2004l . l2006t iTeets et al.ll201ll ). 
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other YSOs that experien ced mass accretion outbursts were found to have similarly enhance d 
or strong X-ray activity (lAudard et al.ll2010t ISkinner et al.ll2006l 120091 ; iGrosso et al.ll2010l ). 
The gravitational potential of YSOs is not deep enough to thermalize plasma to such a 
high temperature; yet the elevated levels of X-ray emission observed from YSOs undergoing 
accretion outbursts (such as V1647 Ori) strongly indicate an intimate connection between 
accretion activity and X-ray emission in these objects. In order to determine where in the 
star/disk system the hard X-rays actually originate, we must ultimately understand the 
origin of the high energy activity during accretion outbursts. To this end, we have reana- 
lyzed all X-ray observations of VI 647 Ori in outburst, in search of temporal evidence that 
might reveal the site(s) and, perhaps, mechanism(s) responsible for its enhanced high-energy 
emission during these events. 



2. Data Sets 



Since the onset of the first outburst in 2004, V1647 Ori has been monitored 17 times 



with three major X-ray ob servatories: 13 time s with Chandra (IWeisskopf et al.l 120021), thre e 
times with XMM-Newton ( Uansen et al.ll200ll ) and once with Suzaku (IMitsuda et al.l 120071 ). 
Six Chandra observations during the first outburst and an XMM-Newton observation during 
the second outburst did not collect enough photons for timing analysis. Table [TJ lists the 
9 data sets used for the timing analysis and the 2 additional data sets added to the folded 
light curve plots (Figs, [fl [2]). Hereafter, individual Chandra, XMM-Newton and Suzaku 
observations are designated CXO, XMM and SUZ, respectively, subscripted with the year, 
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month and day of the observation. 



We reprocessed data with the following calibration versions: DS 7.6 or later for the 
Chandra data, SAS 10.0.0 or later for the XMM-Newton data and version 2.2.11.22 for 
the Suzaku data. The basic analysis fol l ows p rocedures descri bed in earlier pape rs - 



Chandra ( Kastner et al. 



2006 



Teets et al.l 120111 ) . XMM-Newton flGrosso et all 120051 ). and 



Suzaku ( IHamaguchi et al.l 120101 ) . A notable departure from those procedures is that we 
converted event arrival times to the solar barycenteric time system, which differs by up to 
300 sec from the terrestrial time. We generated background subtracted light curves for pho- 
tons with energies in the range of 1—8 keV, binned into 2 ksec (= AT) intervals, using the 
HEAsoftQ analysis package. We developed Python codes for the cross-correlation, chi-square 
and physical model studies. 



3. Period Search 

We identified strong similarities in eleven separate X-ray light curve observations of 
V1647 Ori obtained during the two outbursts with Chandra, XMM-Newton and Suzaku. 
These flux variations — an order of magnitude on timescales of hours — superficially resemble 
stellar magnetic flares, but there is reason to doubt this interpretation, given the spectral vari 



at ion and frequency of flux rises (IGrosso et al.l 120051 ; iKastner et al.l 120061 ; IHamaguchi et al. 



2010|). Fig. [TJ displays all 11 light curves with time offset and flux normalization according 
to the x 2 study described below. We first focus on the light curve with the longest duration, 
obtained with XMM-Newton in 2005 (Fig. [TJ top; the numbering scheme in the following sen- 
tence corresponds with the labels at the top of the figure). The X-ray flux i) stays constant 
for ~20 ksec, ii) rises by a factor of ~5 in ~15 ksec, Hi) keeps an elevated level for 30 ksec 
with marginal spikes and dips, and iv) falls gradually to the original flux level on a similar 
time scale. The XMM-Newton light curve in 2004 apparently matches with — and the 
Suzaku light curve in 2008 starts from iv) and connects to i) — iii). Although the Chandra 
light curves (Fig. [TJ bottom) are of more limited durations and have lower photon statistics, 
each of those light curves also matches with one or more parts of profile z)— iv). 

A cross-correlation analysis provides quantitative support to these qualitative similari- 
ties (see Appendix [A] for details). The XMM-Newton light curves in 2004 and 2005 correlate 
strongly (0.92) when they align at their observation starts and the former is shifted back- 
ward by 6 ksec. The XMM-Newton light curve in 2005 and the Suzaku light curve in 2008 
correlate strongly (0.82) when they are folded by a period of 86 ksec and the Suzaku light 



http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 



-5 - 



curve shifts backward by 30 ksec. 

Given the similarities of these light curves, we search for a more accurate period and set 
of phase shifts that match both the shapes and the timings of all the available light curves (see 
Appendix [B] for details). First, we generate a template light curve with a 1.23 day span by 
combining the XMM-Newton and Suzaku light curves that are shifted in time and normalized 
in flux based on the cross- correlation study. We repeat the template light curve with a 
frequency below 1.45 day -1 , giving a phase offset between 0.0—1.0 and a flux normalization 



for each light curve to account for the long-term variation ( iTeets et al.ll201ll ). Based on our 
preliminary analysis, we also introduce a phase gap between the first and second outbursts. 
We fit the template to all light curves with the least \ 2 method and derive the minimum x 2 
value of 317.9 [x 2 /d.o.f = 1.88 (d.o.f. = 169)] at /„ = 0.98929 day -1 (P =87.3 ksec) and 
&gap — —0.383. The best-fit period is close to the rotation period of the Earth; however, 
because none of the (space-based) observatories whose data are analyzed here obtain data 
on a daily cadence, this period cannot be an artifact of our observing protocol. We therefore 
conclude that VI 647 Ori displays periodic variation of its X-ray emission with a period of 
~1 day. 

One obvious potential origin for the X-rays would be the rotation into and out of our 
line of sight of a localized region of X-ray plasma on the central star. Rises and falls in 
the light curves would correspond to appearances and disappearances of the localized X-ray 
bright spot. The flux transitions take ~20% of one cycle (Fig. [31 Appendix [U]), suggesting 
that the spot has a significant size or height compared with the radius of the central star. 

We assume a uniform circular spot with a tip-cut cone shape and simulate an X-ray light 
curve for each combination of the spot radius, height, latitude and stellar inclination angle. 
We find that no single spot can produce both the low flux interval (0 ~0. 0—0.2, 0.8—1.0) and 
the high flux interval (0 ~0.4— 0.6), so we adopt a bipolar geometry and add a fainter spot 
at the opposite longitude and latitude of the brighter spot whose shape is identical to that 
of the bright spot (see Appendix ICl for details). The best-fit result (x 2 =98.2, d.o.f. =36) is 
obtained with a bright spot to faint spot brightness ratio of 5, with the spots having radii 
of 0.32 and heights of 0.01 i?*, and with the bright spot found at a stellar latitude of 
~ —49°. The stellar inclination — the tilt of the polar axis toward our line of sight — is 
~68°. This model, shown in Fig. [3j reproduces the low and high flux levels and the rise and 
fall of the represented light curve well. The possible dip at $ = 0.5 may be reproduced by a 
partial eclipse of the bright spot by its own accretion flow or the disappearance of the faint 
spot behind the central star if it has a smaller size and higher latitude than those assumed 
in the fit above. Excesses at $ = 0.4 and 0.8 after the flux rise and fall may represent 
asymmetries of the spot shapes or the presence of additional smaller spots. The best-fit for 
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the inclination angle of t he rotation axis is similar to the angle (~61°) estimated from an 



infrared light echo study (lAcosta-Pulido et al.ll2007l ). 



The one-cycle light curve can also be reproduced given plasma configurations that differ 
somewhat from the geometry described in Fig. 3. For example, the faint hot spot — which is 
required in the preceding model so as to reproduce the low flux levels during phases intervals 
$ = 0.0 — 0.2,0.8 — 1.0 — may shift toward the latitude and longitude directions with 
respect to the bright hot spot, or can be replaced by a constantly visible plasma in the form, 
e.g., of a halo around the star. Furthermore, the bright hot spot could really be a complex 
of multiple, smaller hot spots, instead of a uniform single spot. However, in any plasma 
configuration, the flux increase during the phase interval $ = 0.2 — 0.8 constrains the bright 
hot spot (or the envelope of multiple hot spots) to have the approximate size, height and 
latitude described above. 



4. Discussion 



The confidence ranges of the stellar inclination and the latitude of the bright spot (Fig. H] 
left) exclude solutions involving a bright spot in the hemisphere facing us, a bright spot at 
a high latitude, and a pole-on view of the system. This result is consiste nt with modeling 



of strong fluorescent iro n lines observed in Suzaku and Chandra spectra ( IHamaguchi et al. 



20101 ; iTeets et al.l 1201 ll ). which suggests that a significant fraction of hard X-ray-emitting 
plasma is occulted from view. The spot radius can be as large as ~0.5-R*, while its height 
should be lower than ~0.1i?* (Fig. H] right). The spot — shaped like a thin, extensive plate 
- is similar to shape to the mass accretion footpoints of neutron stars, white dwarfs, and 
the Earth's aurorae. 

No periodic variation such as that seen in the X-ray regime has been identified in 
optical or infrared observations of V1647 Ori, most likely because the optical and infrared 
emission mostly comes from the disk. Based on a pre-outburst bolometric luminosity and 
stellar effective temperature, lAspin et al.l (120081 ) roughly estimated the mass and radius of 
V1647 Ori at ~0.8M Q and ~5_R Q , respectively. The stellar radius is slightly larger than the 
distance at which an orbiting body around a O.8M star would have an orbital period of 
about one day@ Thus, the central star must be rotating at a speed close to break-up velocity 
(i.e., rotating at the Keplerian velocity at the stellar radius). 



2 Hcrcaftcr, wc assume a stellar radius of 4 solar radii — the approximate maximum radius that a O.SMq 
star with rotation period ~1 day can maintain without breaking up. 
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Fig. H] right also plots the product of the plasma density squared {n 2 ), the plasma 
filling factor (77), and the cube of the stellar radius (R*), using the plasma emission measure 
determined during the Suzaku observation in 2008. There is no solution at log r]n 2 {R^I '4R & ) 3 
<2x21 cm" 6 . Since 77 5^1 and R* < 4i? , we estimate n >5xl0 10 cm" 3 . Being only a 
lower limit, t his result for den sity — which is similar to those of the densest active stellar 



coronae (e.g., iNess et al.ll2004j ) — may place V1647 Ori in the density regime inferred for 
the footpoints of free- fall accretion (as me asured for isolated T Tauri stars via line ra tios of 
helium- like ions; e.g., Kastner et al.l 120021 for TW Hya; see also iPorquet et al.ll2010l ). The 
magnetic field B should be stronger than B >100 Gauss at the footpoints, in order to confine 
such a dense hot plasma (i.e., if the magnetic pressure is stronger than the plasma pressure; 
plasma (3 = nkT/(B 2 /8ir) < 1). 

The phase gap that we infer between the first and second outbursts may be caused 
by drift of the magnetic poles on the stellar surface or a change in the stellar rotational 
frequency. For the latter case, we can replace the phase gap in our model ephemeris with a 
frequency derivative. In doing so, we find a similarly good solution at a similar frequency 
with a small derivative (see Appendix IB]) . 

The observed X-ray variation of V1647 Ori can be naturally explained by rotational 
modulation of X-ray bright spots. Co ronal active regions can p roduce such rotational X- 



ray modulation (Flaccomio et al.ll2005f ). However, earlier studies (IKastner et al.ll2004l 12006 



Teets et al.ll201l[ ) indicate that large increases in X-ray flux from V1647 Ori during the out- 
bursts are very closely correlated with the (accretion-driven) optical/infrared flux; therefore, 
these X-ray eruptions are most likely accretion-driven as well. The geometrical model de- 
scribed in $3] (Fig. [3]) supports such a model, in that it indicates that the X-ray-emitting 
plasma lies at or very near the foot-points of mass accretion streams at the stellar surface. 

Since the profile of the X-ray light curves did not change remarkably between the ob- 
servations, the intrinsic X-ray luminosity of the hot spots, that is, the accretion-induced 
magnetic reconnection activity, varies on a timescale of a week or longer. Given the large 
overall variation in the amplitude of the X-ray flux from outburst to outburst (and even 
during outbursts) (see Fig. [2]), it is evident that the large-scale magnetic field configuration 
of the V1647 Ori star/disk system is preserved, even as the protostar undergoes dramatic 
changes in accretion rate. We suggest two possible mechanisms that may generate this con- 
dition: (i) the mass ac cretion flow stably disrupts the magnetic fields at the footpoint (e.g., 
Brickhouse et al.ll2010l ). or (ii) the rotational shear between the star and the disk continu- 



ously twists the stellar bipolar magnetic fields (e.g., iGoodson et al.lll997l . see also Fig. E]). 



More evolved YSOs also have faint, hard X-ray emission from hot plasma. This emission 
has usually been explained as due to a blend of emission from multiple micro-flares (e.g., 
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Caramazza et al.l 120071 ). Our result demonstrates that the mass accretion activity also can 



generate hot plasma constantly by sustained magnetic reconnection. The same mechanism 
may also operate on those YSOs with weaker mass accretion activity. 



5. Conclusion 



We have discovered rotational modulation of X-ray emission from the Class I protostar 
V1647 Ori via analysis of 11 X-ray light curves obtained with the Chandra, XMM-Newton 
and Suzaku observatories during this YSO's two recent mass accretion outbursts. Based on a 
cross-correlation study and period search, we determined a rotational period of ~1 day, with 
either a phase gap apparent between the two outbursts or frequency variation. The single- 
cycle light curve can be successfully reproduced by emission from two hot spots on opposite 
poles on the stellar surface. The star rotates with a period of ~1 day, close to the break-up 
velocity for a 0.8 M & star with a radius of ~4 R & . The hot spots likely cover significant 
fractions of the stellar surface, while the height (~0.01 R*) may be negligible compared with 
the stellar radius. The hot spot size and the plasma emission measure indicate relatively 
high plasma density (>5xl0 10 cm -3 ), also pointing to an origin in accretion hot spots. This 
result clearly demonstrates that hard X-ray-emitting plasma can be present in long-lived 
accretion footprints at the surfaces of protostars, and thereby constrains the geometry of 
magnetospheric accretion in early (Class I) protostellar evolutionary stages. 



This work is performed while K.H. was supported by the NASA's Astrobiology Institute 
(RTOP 344-53-51) to the Goddard Center for Astrobiology (Michael J. Mumma, P. I.). J.K.'s 
research on X-rays from erupting YSOs is supported by NASA/GSFC XMM-Newton Guest 
Observer grant NNX09AC11G to RIT. This research has made use of data obtained from 
the High Energy Astrophysics Science Archive Research Center (HEASARC), provided by 
NASA's Goddard Space Flight Center. 
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A. Cross-Correlations of the Light Curves 



We cross-correlate light curves of XMM040404, XMM050324, and SUZosioos, which have 
good photon statistics. We define the cross-correlation index r, 

E[(^ -x)(yi-d-y)] 



where X{ and yi are count rates of the 2-th bins of two light curves, d is a delay in units of 
2 ksec time bins, x and y are averages of bins that contributes to the cross-correlation. The 
index r ranges between —1 and 1. Two light curves are identical when r — 1. 

We first cross-correlate XMM040404 (xi) with XMM 0503 24 (yi) (Fig. [6] left). There is a 
strong correlation of r = 0.92 at d = 3 (6 ksec). 

Since SUZ 8ioo8 and XMM050324 have similar observing durations, the number of bins 
that contributes to the cross-correlation becomes smaller at larger delays. The r index 
fluctuates strongly in these regions. We thus require the number of contributing bins to 
be at least 10. When SUZ 081008 (xj) is shifted backward against XMM 050 324 (yi), there is 
a strong correlation of r = 0.91 at d = 28 (56 ksec). When XMM 050 324 (xi) is shifted 
backward against SUZosioos (yi), there is a strong correlation of r = 0.81 at d — 15 (30 ksec). 
This result indicates that these two light curves folded at a certain period also match well. 
We thus fold both light curves by 40—46 bins (80—92 ksec), average overlapped bins with 
weighted mean values, and cross-correlate them. We found that the r index is a maximum 
of 0.82 when these light curves are folded by 43 bins (86 ksec) and XMM050324 (xi) is shifted 
backward by d = 30 (60 ksec) against SUZqsioos (yi) (Fig. M right). 



The cross-correlation tests the similarity of two light curves, but it does not consider 
the time interval between them. We therefore search for an ephemeris that satisfies both 
shapes and timings of all the X-ray light curves, including those detected with Chandra. 

We first generate a template light curve from the XMM040404, XMM050324 and SUZosioos 
light curves. Based on the cross-correlation study in Appendix |AJ we shift XMM040404 and 
XMM 050 324 backward by 18 and 15 bins (36 and 30 ksec) against SUZ 08100 8, respectively. 
We normalize these light curves based on the average count rates of overlapped bins and 
average them at weighted mean values. After artificially increasing the time resolution of 
the averaged light curve by 100 times with linear interpolation, we smooth it with a Gaussian 



r = 




(Al) 



B. Search for the Best Ephemeris with the \ 2 Test 



- 10 - 



function with a = 2 ksec so as to minimize statistical noise. The resulting template light 
curve [L temp (ti)] spans a time period of 1.23 day. 

We define the ephemeris as: 

$(T) = % + fo (T - T ) + % ap H(T - T gap ) (Bl) 

where f is frequency, T is the time origin, $ is the phase at T = T , § gap is the phase gap, 
H is the unit step function and T gap is the time of the phase gap. We then assign phases, 
*&(ti), to the template light curve according to 

*(fi) = foU (B2) 

and discard bins at \&(tj) > 1. 

We sum bins of the template light curve within the corresponding phase range of each 
bin of a measured light curve, that is, 

n 

Itemp(Ti) = ^ L temp (tj ) (B3) 



where m and n satisfy 

$'(T,) = $(Ti) - L<f(T)J (B4) 
A$ 

*(t m ) ^ $ , (T j )-y<$(Wl) (B5) 
A$ 

tf(i B ) ^ $'(T,) + — < tf(i B+1 ) (B6) 

where [^J is the floor function and A$ = /oAT. We then normalize Itemp{Ti) for each 
measured light curve by the normalization factor n G & s that gives a minimum x 2 value, and 
derive a \ 2 value for each set of $ , "^Vp, and f . 

XI L bs(Ti) Itemp(Ti) I AL D bs(Tj) 2 /T37\ 
2 / Lobs(Ti) — n o b s Itemp(Ti) 2 /DO > 

x = ^ ( al^ } (B8) 



obs 



We search for the minimum x 2 value in the ranges of ^ $o> $gap < 1 and /o ^ 0.82 day 1 . 
We find a minimum x 2 of 317.93 \x 2 /d.o.f. = 1.88 (d.o./. =169)] at /„ = 0.98929 day" 1 
(P =87.3 ksec) and § gap = —0.383. The best ephemeris is hence expressed as 

$(T) = 0.98929 (T - 13453.09075) - 0.383 H(T - 14300) (B9) 
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where T is the truncated Julian date. In this formula, we re-define the $ origin at an 
intermediate point between the fall and rise timings of the single cycle light curve (see 
Appendix IC]) . 

The phase gap is empirical. The phase change may instead be explained by a frequency 
variation. We hence modify equation f ]Blj) to allow for a frequency derivative. 

$(T) = $ + /o (T-T ) + /i (T-T ) 2 /2 (BIO) 

We find a similarly good fit of \ 2 = 317.22 at f =0.98932 day" 1 and fx = -1.89xl0~ 6 day" 2 
at T = 14747.6758097 day. Figs. [7] and M show light curves folded with this ephemeris. 

We generate a single cycle light curve from XMM 40404, XMM050324 and SUZqsioos- We 
fold these light curves with equation (IB9j) and bin all light curves with AT =1984.9 second, 
so that one whole light curve requires exactly 44 bins. We normalize and average these light 
curves in a manner identical to that which produced the template light curve. 



C. Physical Model 



We consider an X-ray point source sitting at longitude 8, latitude and height h from 
the stellar surface, viewed from the rotation axis at inclination angle i (see Fig. [9] left). We 
define the 9 origin as the longitude that crosses the opposite side of the central star from 
the Earth at $ = 0. The point source appears in view at a longitudinal difference A6 e from 
this (#=0) reference point (Figs. M right) that satisfies: 



cosA6 P 



cosi sincf) + yjl — 



R,+h< 



sini coscf) 

The normalized X-ray light curve of the point source is: 

Tpoint^) = (0^~<$ e ) 

= 1 ($ e ^ E ^ 1 - $ £ 
= (1-$ K <E<1) 



(CI) 



(C2) 
(C3) 
(C4) 



where 



2tt 



[$ + £J and $ e = 2nA6 e 



X-ray spectra of V1647 Ori indicate that the X-ray plasma is opt i cally thin fjKastner et al 



20061 : iGrosso et al.ll2005l : lHamaguchi et al.ll2010l : lKastner et al.ll200i fleets et al.ll201lh . This 
means that any portion of the X-ray plasma with a finite size can be seen once the portion 
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emerges from the stellar rim, such that the observed flux at phase $ is the integrated emission 
from the plasma emerging above the stellar rim, i.e., 



= J J J ^ P omt($, 9, 0, h, i)S(6, (j), h)(R* + hycos<pd9d<pdh (C5) 

where S(9, (ft, h) is the X-ray source distribution of the plasma. 

For simplicity, we assume a uniform plasma with a conical shape with angular radius 
and height h', standing upside down on the stellar surface, such that the base position 
of the cone is "exposed" (Fig. [TO]) . It is observed at unit flux when in view, such that the 
source distribution is described as 

S(9,(f),h) = 1/V (inside) (C6) 
= (outside) (C7) 

where V is the plasma volume. We define the longitude and the latitude of the cone axis as 
9' and <j)', respectively. A narrow latitudinal strip at + A<p' (|A0'| < r*) has a half width 
A9 S that satisfies 

. cos r* - sin(<f)' + A0')sm0' 

cosAu s = — (L8) 

cos(<p' + A(j)')cos(J)' 

We numerically integrate equation (105|) for this plasma, i.e., 

JO J<j>'-r t Jd'-Ae s * 

The low and high flux phases in the single cycle light curve cannot be reproduced by 
any single spot. We therefore assume two spots with identical shapes in a dipole geometry, 
sitting on opposite sides of the central star. The X-ray light curve is, then, expressed as 

F = g b F conc (<5>, 9', <f)',n,ti,i)+g f F conc (<5>, 9' + tt, -<f)', r*,ti,i) (CIO) 

where gj, and gj are un-occulted fluxes of the bright and faint spots, respectively. The x 2 
value is a minimum (98.2) when <p' = -49°, r* = 18°, h! =0.01i2*, i =68° and g f /g b = 0.20. 
In this fit, we fixed the longitude 9' at 0, considering the definition of the $ origin. Fig. [3] 
illustrates this best-fit model. 

The plasma emitting volume is, 

Ve = VV (Cll) 

= r] / (R* + h) 2 ttdh (C12) 



2vrr/(l - cos r*)ti 2 



[3Rl + 3R*ti + h' 2 } (C13) 
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where rj and f2 are the plasma filling factor and the solid angle of the cone, respectively. To 
constrain the plasma density, we combine the standard relation for plasma emission measure, 
E.M. = n 2 V e (n: plasma density), with equation fl013j) . 

2/ R* \3 E.M. (n-\A\ 

VU [ 4R^ } ~ 2vr(4R ) 3 (l - cos n)h*{l + K + h 2 j3) [ ' } 

where h* = h'/R*. We calculate the right side of the equation for each combination of r 
and h * and E.M. = 1.9 x 10 54 cm -3 during the Suzaku observation in 2008 ( jHamaguchi et al. 



20101 ). Fig. H] right shows contours of values for this parameter. 
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Table 1. Analyzed Data Sets 



Abbreviation 


Observation ID 


Start Date 


Exposure 


Duration 


Net Count 








(ksec) 


(ksec) 


(counts) 


First Outburst: 












CXO040307 


5307 


2004 Mar. 7 


5.5 


5.6 


60 


CXOo40322t 


5308 


2004 Mar. 22 


4.9 


5.0 


10 


XMM 40404 


0164560201 


2004 Apr. 3 


29.1 


37.0 


1321 


XMM 50324 


0301600101 


2005 Mar. 24 


79.2 


89.7 


1557 


CXOo50411 


5382 


2005 Apr. 11 


18.2 


18.4 


85 


Second Outburst: 












CXO080918 


9915 


2008 Sept. 18 


19.9 


20.2 


455 


SUZ0810O8 


903005010 


2008 Oct. 8 


40.4 


77.5 


1275 


CXO08I127 


10763, 8585 


2008 Nov. 27 


20.0, 28.5 


77.5 


197,143 


CXO090123 


9916 


2009 Jan. 23 


18.4 


18.6 


240 


CXO090421 


9917 


2009 Apr. 21 


29.8 


30.2 


258 


XMMl00228t 


0601960201 


2010 Feb. 28 


33.5 


34.0 


163 



Note. — Abbreviation: CXO — Chandra, XMM — XMM-Newton, SUZ — Suzaku. f The 
datascts are not used for the timing analysis. Net count: Background subtracted photon 
counts between 1—8 keV. Photon counts of all the available instruments are summed for 
XMM-Newton and Suzaku. 
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Fig. 1. — Light curves of V1647 Ori between 1—8 keV obtained with the XMM-Newton, 
Suzaku (top) and Chandra (bottom) observatories. These light curves are folded and normal- 
ized according to the best fit ephemeris (fo =0.98929 day" 1 , $ ffop = —0.383; see Appendix IB]). 
The upper and lower horizontal axes in each panel show the rotational phase and time in 
second from the phase origin, respectively. The phase origin is set at the middle of the flux 
fall and rise. Light curves in thin colors are repeats of thick ones in the same colors. Points 
with dotted error bars are not used for the best ephemeris search. Solid purple and green 
lines in the top and bottom panels are the template light curves for the \ 2 analysis. Curves 
with narrow widths indicate phase intervals that have been disregarded, as per the algorithm 
described in Appendix [Bj The arrows at the top depicts the phases of variation mentioned 
in Section [3j 
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Fig. 2. — Same as Fig. HJ but with axes in energy flux units. 
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Phase O 

Fig. 3. — Single cycle light curve created from the XMM-Newton and Suzaku light curves 
overlaid with the best-fit result of the symmetric bipolar spot model [red: total, green: bright 
spot, blue: faint spot). The pictures on the top depict locations of hot spots at corresponding 
phases. See Appendix O for the numbers and letters in the label. The subscripts b and f 
stand for the bright and faint spots, respectively. 
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Fig. 4. — Confidence ranges of parameter values obtained from model fitting (Fig. [3]) - 
left: the stellar viewing inclination vs. the latitude of the bright spot; right the radius 
vs. height of the hot spots, in stellar radii. In both panels, the black, red and green 
contours show confidence levels of 66% {A X 2 =2.3), 90% (A^ 2 =4.6) and 99% (Ax 2 =9.2), 
respectively. The blue stars show the best-fit result. The right panel also shows contours 
of log r]n 2 (R* /4:Rp ) ) 3 = 21, 22, 23, a nd 24 in solid orange lines as obtained from modeling 
Suzaku data (IHamaguchi et al.ll2010l see Appendix [C] for the derivation.) 
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Fig. 5. — A possible geometry of the V1647 Ori system (mechanism ii). Differential rotation 
between the star and the disk shears the stellar bipolar magnetic fields, and the magnetic 
fields twist and continuously reconnect. Matter, accelerated to v ^2000 km s -1 by the 
magnetic reconnection, collides with the stellar surface, th ermalizes to kT ~4 keV ( T ~50 
MK) and emits hard X-rays (See e.g., the equation (9) of ICalvet &: Gullbringlll998l ). The 
opposite magnetic pole may have lower mass accretion and emit weaker X-rays. 
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Fig. 6. — Cross-correlation of XMM050324 with XMM040404 (left) and SUZ 8ioos (right) when 
these light curves are folded by 43 bins (86 ksec). 
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Fig. 7. — Normalized light curves (top: XMM-Newton, Suzaku, bottom: Chandra) between 
1—8 keV folded with the ephemeris in the equation ( IBlOj) . which assumes a frequency deriva- 
tive instead of a phase gap. See the caption of Fig. CD for the other details. 
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Fig. 8. — Same as Fig. but with axes in energy flux units. 
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Fig. 9. — Definitions of the coordinate system, 6, <p, h and i (left) and of A8 e when h = 
(right top: edge-on view, right bottom: pole-on view). 




Fig. 10. — Definition of the parameters 9', <fr', h', A<p', Ah and A9 S for the conically shaped 
plasma model. 



